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Abstract. We investigate the distribution properties of the fractional Levy motion. We 
consider separately the cases Q < H < 1/2 {short memory) and 1/2 < iJ < 1 {long memory), 
where H is the Hurst parameter, and present the asymptotic behaviour of the distribution 
density of the process. Some examples are provided, in which it is shown that the behaviour 
of the density in the cases < iJ < 1/2 and l/2<iJ<lis completely different. 

Keywords: Fractional Levy motion, distribution density, saddle point method, sub-exponential 
distribution. 
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1. Introduction 

In this paper we consider the distribution properties of the fractional Levy motion (FLM, in 
the sequel). Various versions of the FLM have been used in a number of recent publications in 
order to interpret some experimental data. Apart from the rigorous mathematical definition, 
some modifications of the FLM are derived from the physical point of view, see, example, 
|H99) . |CSC09| . |ES12] . The FLM driven by an a-stable Levy process is used as a model for 
describing sub-diffusive effects in physic and biology (see |BW10| . |WBMW05| ). signal and 
traffic m odeling |LLHDn2|, [KM96| . jDMQS], finance |BSVn8) . geophysics |PP94| . [WoH] . 
|GLS06) . We refer to |W05) for the discussion in which type of problems the FLM gives 
an adequate description for the observed phenomena. In the papers quoted above it was 
shown that the respective phenomena existing in nature can be better described by models, 
containing the FLM rather than the fractional Brownian motion (FBM). Finally, we refer to 
|CGOO| for simulations of the FLM, which can be convenient in practical problems. 

Similarly to the FBM, the FLM can be defined in two different ways: via the Mandelbrot- 
Van Ness representation (see |BCI04] and |M06] ). or via the Molchanov-Golosov represen- 
tation (see |BM09) ). We also refer to |ST94] for a bit different definition of the fractional 
stable motion. These two representations, being equivalent in the Gaussian setting, in the 
Levy setting lead, in general, to different processes; see |TMllj . Note that, in contrast to 
the FBM, in some cases the FLM can even be a semi-martingale ( |BP09) . |BLS12] ). 

In this paper we focus on the FLM defined by the Mandelbrot- Van Ness representation, 
i.e. 



f{t, s)dZs 



where Z^, s G R, is a (two-sided) real- valued Levy process, H & (0, 1) is the Hurst parameter, 
and 



f{t,s) 



\H~l/2 / \H~l/2 



t,s eR, 



[1.2) 



T{H + 1/2) I' 

where x+ := max(a:, 0). This definition gives a particularly important representative of the 
class of so-called moving- average fractional Levy motions. Since the FLM, according to the 
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survey above, is an adequate model to some phenomena in nature, it would be appropriate 
to investigate deeply its properties. In particular, knowledge of the distribution properties 
of the FLM would naturally make it possible to solve various problems related to statistical 
inference, simulation, etc. 

In this paper we concentrate on the asymptotic behaviour of the distribution density of 
the FLM. In contrast to the FBM case, the study of the distribution density of the FLM 
is much more complicated. In the recent paper |KK11| we presented the investigation of 
the distribution density of a fractional Levy motion in the following cases: (i) H = 1/2, 
which means merely Zf^ = Zt, and (ii) 1/2 < H < 1, which corresponds to the so-called 
long memory case, see Definition 1.1 in |M06| . Both cases can be treated in a unified 
way using a general result about the asymptotic behaviour of distribution densities of Levy 
driven stochastic integrals with deterministic kernels, see Theorem 2.1 in |KK11| . Since this 



theorem requires the respective kernel to be bounded, and the kernel (II. 2p is unbounded 
when < H < 1/2, the case (iii) < H < 1/2 cannot be treated in the same way as in 
|KK11| . and thus requires a completely different approach, which we present below. Further, 



we show that there is a substantial difference in the behaviour of the density in the cases 
(i), (ii) on one hand, and the case (iii) on the other hand, namely, the distribution density 
in these situations exhibits absolutely different types of the asymptotic behaviour. We also 
emphasize, that in contrast to the situation studied in |KK11] . for the case 0<if<l/2we 
do not require the existence of exponential moments of the tails of the Levy measure. 

Let us outline the rest of the paper. In Section [2] we set the notation and formulate our 
main results. Theorems 12.11 and 12.21 In Section [2] we formulate the general result. Theo- 
rem I3.lt on the asymptotic behaviour of the distribution density of Levy functionals, on 
which Theorems 12.11 and 12.21 are based on. In Section H] we give two examples which illus- 
trate the effects that may happen in the "extremely heavy-tailed" case, i.e. when condition 
(I2.14P (see below) fails. In Appendix we give some supplementary statements: the necessary 
and sufficient condition for the integral (11. II) to be well defined, and the condition for the 
respective distribution to possess a density. 



2. Settings and the main result 

Let Zt, t G R, be a real-valued (two-sided) Levy process with the characteristic exponent ijj, 
which means that Z has stationary independent increments, and the characteristic function 
of an increment is given by 

The characteristic exponent ijj admits the Levy-Khinchin representation 

ip{z) = iaz - bz'^ + / (e™^ - 1 - 22m1{|„|<i}) /i(rfM), (2.2) 



where a G M, 6 > 0, and yu(-) is a Levy measure, i.e. J^{1 A u'^)fi{du) < oo. To exclude the 
trivial cases, we assume that 6 = and /i(El) > 0; that is, Z does not contain a diffusion 
part, and contains a non-trivial jump part. To simplify the notation we also assume without 
loss of generality that Zq = and a = 0. 

The integral (11. ip is defined as a limit in probability of the respective integral sums, see 
|RR89) . Section 2. When H ^ 1/2, the necessary and sufficient condition for this integral to 
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be well defined is 

/ |n|2/(3-2^^)^(rf^) < oo, (2.3) 
J|«|>i 

see Proposition lA.ll in Appendix A. Furthermore, it will be shown in Proposition IB.ll (see 
Appendix B) that under the same conditions and our standing assumption 

/x(R) > 0, (2.4) 

the integral fll.ip possesses for any t ^ a distribution density, which we denote by Pt{x), 
and, moreover, pt G C^(R). Note that in the Levy case, i.e., for H = 1/2, available sufficient 
conditions for the existence of the density pt G C^(R) are much stronger, see, for example, 
[HW42] for the Hartman- Wintner condition. 

An important feature of the process Zf^ is that one can explicitly write its characteristic 
function (t){t,z) := Ee'^^? (cf. Theorem 2.7 in |RR89| ): 

0(t,z) = e*(*'-^), (2.5) 

where 



\l/(t, z) := exp 




^^-izf{t,s)u _ ^ ^ .^^^^^ s)m1|„|<i) li{du)ds 



2 G R, t > 0. (2.6) 



Observe, that if the measure /i possesses exponential moments, the function $(t, z) can be 
extended with respect to z to the complex plane. Moreover, one can see (cf. Section 13. 3p 
under the assumptions that yu(R+) > 0, the function 

H{x,z) := izx + '^{l,z) 

has a unique critical point ^{x) on the line iR. Put 

D{x) = H{x,tax)), K{x) = g^H{x,tO (2.7) 



and 



Mfc(0 := / u'=e«"/i(rfu), k>2. 
Jr 



Fix to > 0. 

Theorem 2.1 ([KKll]). Let Zj^ , 1/2 < H < 1, t > to, be a FLM defined by where 
Zf is a Levy process with the associate Levy measure fi. Suppose that the conditions below 
hold true: 

1) /i(R+) > 0; 

e^yfx{dy), forallxeR] (2.8) 



'\y\>i 

3) 37 G (0, 1) such that M^{^) < M|(70 as ^ ^ 00; 

4) In V 1) + lnlnM2(0 « ^ ^ +00. 



Then the transition probability density pt{x) of Z^ exists, pt G C'^ , and satisfies the 
asymptotic relation 

Ptix) ^ ^tDjxt-^-^/^) (t,x) e [to, 00) X R+, (2.9) 

^ ' ^2Tit^^iK{xt-H~i/2) > V , ; L 0, ; +, v ; 

where D and K are given by (12. 7p . 
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To formulate the result for < < 1/2 we need a bit more notation. Let f{s) := /(I, s), 
see (11. 2p for the definition of f{t, s). Observe that f{s) is monotone decreasing on (— oo, 0), 
monotone increasing on (0,1], and the range of / restricted, respectively, to (— oo,0) and 
(0, 1], is (— oo, 0) and [lyj^iijfj, +oo). In addition, the derivative /'(s) is well defined and is 
continuous on (— oo,0) and [0, 1). Hence, we can put 



ye (-00,0) U 



^ '-00 



£(y):=^^-' ~ nr^{y))' ^'"^ ^ [r(i/+i/2)' ' ^7' ^2.10) 
[O, otherwise. 

Note that i{y) is non-negative if ?/ > 0, and is negative otherwise. Define 

m(r) := / -i (-] fi{dy), r > 0. (2.11) 

J-ooV \yj 

Recall that (see Definition 4 in |K189) ) a function g : R_|_ — )■ 1R_|_ belongs to the class Ld of 
sub- exponential densities, if g{x) > for large enough positive x, and 

lim ~ = 1 for any y G R. (2.12) 

x^+oo g[x) 

Fix t > 0. 

Theorem 2.2. Let , < H < 1/2, be a FLM defined by (11. ip . where Zf is a Levy process 
with the associate Levy measure fi. Suppose that 

/i(R_) > 0. (2.13) 

Then there exists the transition probability density pt{x) of FLM (11. ip . and pt G C^(R). 
Further, 



then 



with 



< oo, (2.14) 

|m|>1 

Pt{x)r^CH[ I |n|2/(i-2^)^(dn)^x~(3-2//)/(i-2i/)^ a;^+oo, (2.15) 



2 / / \\-2/(l-2-ff) 

CH = -, — r iJ+1/2)) . (2.16) 



1 - 2H 

a) if fl2.14p fails, but (12. 3 p holds true, and m G Ld, then for every positive t 

Pt{x) ~ t=^/2-/f^(^i/2-//^^^ ^ ^ ^2.17) 

Remark 2.1. a) Apparently, the relation (12.171) holds true when both x, t tend to +oo in 
such a way that xt~^~^^'^ — ?■ +oo. To prove such an extension of Theorem 12. 2^ one should 
have an extension of the Theorem 3.2 in |K189) which applies to a family of random sums 
with the variable distribution of the number of summands. 

b) Clearly, results similar to (I2.17P and (I2.15P can be formulated for x — oo. In that 
case one should assume /i(R+) > instead of yu(R_) > 0. 

To illustrate the crucial difference between the cases treated in Theorem 12.11 and 12. 2[ 
consider two particular examples from |KK11| which concern the case H > 1/2. First, let 
the Levy measure fi of the Levy noise Zt in (II. ip be supported in a bounded set. Then (see 
Corollary 5.1 and Corollary 5.2 in |KK11| ) there exists a constant c*(/i), defined in terms of 
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the Levy measure fi only, such that for any constants Ci > c*(/i) and C2 < c*(/i) there exists 
y{ci, C2) such that for x > y{ci, 02)^^'^^^^ we have 

Mx) ;-"^P(-W7§P^l^(^)y' (2.18) 



< 



exp y r{H+i/%tH-i/2 1^ (4^+1/2)^ 



Similar statement is available as well when the tails of the Levy measure admit the following 
super-exponential estimates: for u large enough 



^e-^"" < f,{[u, +oc)) < Q(n)e-^"', (2.19) 

where 6 > and /3 > 1 are some constants, and Q is some polynomial. In this case, instead 
of f l2.18|) we have for any constants ci > c*(/i) and C2 < c*(/i) 



P,W i'""} ^'""'S'"^' S'Tii (2.20) 

for X > y{ci, C2)t^^^^'^ (again, c*(/i) is defined in terms of the Levy measure /i only). 

Comparing fl2.18p and f l2.20p . we see that the asymptotic behaviour of the tails of the 
Levy measure fi is substantially involved in the estimates for pt{x). The case < H < 
1/2 is completely different. In particular, if fl2.14p holds true, then pt{x) satisfies (I2.15p . 
where the right-hand side is even independent of t, which is an in interesting and quite 
an unexpected fact. We also emphasize, that under (I2.14p the polynomial "shape" of the 
expression in the right-hand side of fl2.15p does not depend on /i, and the only impact of fi is 
represented by the multiplier \u\'^^^^~'^^^ n{du). This means that in the case < H < 1/2 
the asymptotic behaviour of pt{x) "mostly" does not depend on fi. However, when /i is 
"extremely heavy-tailed", there still remains a possibility for the density pt{x) to be more 
sensitive with respect to both the Levy measure fi and the time parameter t. The dichotomy 
between the "regular" case (when f l2.14p holds) and the "extremely heavy-tailed" case (when 
f l2.14p fails) is illustrated in Section |4] below. Such a dichotomy can be informally explained 
by the competition between the impacts of the kernel f{t, s) on one hand, and of the measure 
IJ, on the other hand. 



3. Proofs 

3.1. General theorem. Before we proceed to the proofs, we formulate a central analytical 
result on the behaviour of the inverse Fourier transform of a certain class of functions. This 
result plays the key role in the proofs of Theorems 12.11 and 12.21 

Let / C R be some interval and T be some set of parameters. Consider a function 
/ : T X / ^ R such that 

f{t,s)ds< 00, teT, (3.1) 

and a family of subsets C(t, s) C R, t G T, s G /. For a given Levy measure fi and a function 
/ satisfying (13. ip let 

^{t,z):= [ [ (e-^"^(*'*)"-l + iz/(t,s)nl|«l<i)/i(rfu)cis, zER. (3.2) 
J I Je{t,s) 



Qt 
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Our aim is to investigate the asymptotic behaviour of the function (provided it exists) 

(x) = (27r)-i / e-'""+*(*'")d^, t>0, x E R, (3.3) 
Jn 

as t and x tend to infinity in some appropriate regions. Clearly, when C(t, s) = R, the 
function qt{x) is nothing else but the distribution density of the Levy functional 

Yt = jj{t,s)dZ,, (3.4) 

where Zg is the Levy process associated with measure fi, without a drift and a Gaussian 
component. 

Assume now that 

HI $(t, ■) possesses an analytic extension to C-|_ := {z G C : Imz > 0}. 
Under this assumption the function 

H{t,y,z):=tzy + ^{t,z) (3.5) 
is well-defined for z E C^. Observe, that -^H{t, x,i^) > 0, and 
d 

—H(t,x,0) = 0, and H{t, x, i^) +oo, ^ — > +oo. 
Therefore there exists the solution $,(t,x) to the equation 

-^H{t,x,tO = 0, (3.6) 

which is unique for every a; > 0. 

To formulate the result we need some extra notation: 

T>{t, x) := H{t, X, i^{t, x)), X{t, x) := M2{t, C{t, x)), (3. 

and 



^kit,0-=i^'^it,tO, k>l, (3.7) 



Q{t,z,B):= / / il — cos{f{t, s)zu)] ^{du)ds. 

J I J {u: f{t,s)uGB,u£e{t,s)} ^ ' 

For a subset ^1 C T x R+, define 

T := {t : 3a; G R+, (t, x) E A], S := {(t, : 3(t, x) E A, {t, = (t, ^(t, x))}. 

Finally, suppose that the function 6' : T — t- (0, +oo) is bounded away from zero on T, and 
the function x : T — > (0, +oo) is bounded away from zero on every set {t : 6(t) < c}, c > 0. 

Theorem 3.1. For the function f satisfying flL2p and the Levy measure fi let ^{t,z) be 
defined by (13. 2p . Assume that HI and the conditions below are satisfies: 

H2 M4(t,0 « M2(t,0, e{t)+^^ oo, (t,0 G 03; 

H3 For e{t) +^^oo, {t, {) G 3; 

In V + In ((in ((1 V X~\t))M,{t, 0) ) V l) « In^(t) + xi^ 

H4 There exist R> and 6 > such that 

e(t,2,R+) > (1 + 5)ln(x(t)|2|), tE7, \z\>R; (3.9) 
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H5 There exists r > such that for every e > 0, 

inf e{t,z, [rx(t),+oo)) > c9{t)({ex{t)f A l), t G T, c> 0. 

Then the function qt{x) given by (13. 3p is well-defined, and satisfies 

~ -=4==e^^*'^\ e{t)+x->oc, {t,x)eA. (3.10) 
y 27rX(t, xj 



Up to some straightforward and purely technical modifications, the proof of Theorem 13.11 
coincides with the proof of Theorem 2.1 in |KK11] . and therefore is omitted. Here we only 
remark that the proof is based in the appropriate modification of the saddle point method, 
see |Co65) for details. 

3.2. Outline of the proofs. One can prove Theorem 12.11 using a simplified version of 
Theorem 13.11 with C(t, s) = R, as well as the scaling property of the function f(t, s): 



/(t, .) = t^-^V (^) , (3.11) 



see |KK11| for the detailed proof. 

Let us turn now to the proof of Theorem 12.21 To make the proof of Theorem 12.21 more 
transparent, we first sketch its main idea. In particular, we show how Theorem 13.11 applies 
in the situation when the function f{t, s) is unbounded. 

According to Theorem 2.7 in |RR89] . the characteristic function (f){t,z) (cf. (12. 5p ) can be 
decomposed for any fixed A > as 

where 

Mt, z):= [ [ (e-/(*.^)" - 1 - zzfit, s)nl|„|<i) fi{du)ds, (3.12) 

J-oo J{u:uf(s/t)<X} 

iP2{t,z) := f I {e''^^''''>'' -1) fi{du)ds + iza{t), (3.13) 

^— oo J {u:uf{s/t)>X} 

a{t):=j j f{t,s)ul\u\<Mdu)ds = t^+^/'^a{l). (3.14) 

J~oo J {u:uf{s/t)>\} 

In the last identity we have used the scaling property (13. lip of the kernel /. 

The function |</)i(t, z)| is integrable w.r.t. z for any t > 0, see Remark IB. 2 1 in the Appen- 
dix B below. Then there exists the distribution density 

Ptix):=^ I e-'^>i(t,^)c?^, (3.15) 



27r jR 

and thus the required density Pt{x) can be written as the convolution 

Pt{x) = {Pt*Pt){x), (3.16) 

where Pt{dy) is the probability measure corresponding to the characteristic function 02 (t, z). 
Define the measure Mt{dy) by the relation 

f g{y)Mt{dy) = [ [ g{f{t, s)u) ^i{du)ds, 

Jm. J -oo J {u:uf(s/t)>\} 
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where g is an arbitrary measurable and bounded function. Then, up to the shift by a(t), the 
measure Pt{dy) is equal to the distribution of the compound Poisson random variable with 
the intensity of the Poisson part equal to Mt{dy). In other words, 

Ptidy) = 5ait){dy) * I e-^^*W<5o(ciy) + e"^^'(^) l^Mfidy)] , (3.17) 



k=l 

r*k/ 



where (dy) is the /c-fold convolution of Mt{dy). 

It follows from the scaling property (13. lip that Mt(R,) = tA with 

A= f I fi{du). (3.18) 

^-oo J {u:uf{s)>X} 

Furthermore, we prove that the measure Mt{dy) is absolutely continuous with respect to the 
Lebesgue measure with the density 

mt{x) = t'/^-''mit'/^-''x)l^,^^,H-i/2}, (3.19) 

where the function m is defined by (12. lip . Then (I3.16P can be written in the form 

Pt{x) = e-'^'Ptix - a{t)) + / pt{x - y)pt{y - a{t))dy, (3.20) 

where 



R 



Mx):=e-^^Y.^- (3-21) 



k=l 

Clearly, pt is the density of a random sum with the distribution of one term represented 
by rrit. Suppose that the function m is sub-exponential. Then it follows from Theorem 3.2 
in [K189] that the density pt is sub-exponential as well, and 




where we used that rnt{x) dx = Mt(R) = tA. 

To estimate pt{x) we apply Theorem 13. II Namely, in Proposition 13. II below we show that 
for a given e > there exists y{e) > such that 



exp 



'1 + e)x , X \ _ , , / (1 — e)x , X \ u mo 

In ITTTTT^ < Pt{x) < exp -^-^ In -— - , xr^-'/' > y{e) 



tH-l/2 fH+1/2 I — -fiv-^/ — ^^^t- 1 fH~l/2 fH+1/2 



(3.23 ) 

Since a sub-exponential function decays slower than any exponential function (cf. |K189) ). 
the term rrit dominates both pt{x) and the integral term in (I3.20p . In such a way, (I3.20p . 
(I3.22P and (I3.23P provide the required relation fl2.17p . 

Let us summarize the idea explained above. The distribution of Zf^ is decomposed in two 
parts. For one part, the distribution density is controlled by means of the respective version 
of the saddle point method, while for the other part the distribution can be evaluated in the 
form of the series of convolution powers with the explicitly given law of the first summand. 
Similarly to Theorem 2.1 in |KK11| . Theorem 13.11 provides a flexible version of the saddle 
point method, which is applicable to a wide variety of integrals of the form (13. 3p . Thus one 
can expect that the approach presented above can be extended to other processes of the form 
(II. ip with unbounded kernels f{t, s). To keep the exposition reasonably tight, in this paper 



ASYMPTOTIC BEHAVIOUR OF THE DISTRIBUTION DENSITY OF THE FLM 9 

we do not investigate this possibility in the whole generality, and restrict our considerations 
to the important particular case of the FLM with < < 1/2. 

3.3. Properties of pt{x). 

Proposition 3.1. Under for a given e > there exists y{e) > such that pt{x) 

satisfies (I3.23p . 

Proof. We use Theorem O with e{t) = t, x{t) = t^'K T = [to, oo), / = (-cx),t], e(t, s) := 
(— oo, j^jj]^], and A = [to, oo) x R_|_. Observe that condition HI is automatically satisfied. 
Recall that / possesses the self-similarity property (13. lip . Then 

Mk{t,C) = x'{t)tMk{x{t)C), (3.24) 

where 

M.(0 ■= ^H(l X <(1 X)) = l^-r ^^'"^ " ^ = 

JYL,(Cj. g^,M[i,x,^C[i,x)) <^jX^^k^.cN{dv), k>2. 

Here N(dv) := J^^fls{dv)ds, and jlsidv) is the image measure of fi(du) under the mapping 
u (— !■ f{s)u. It can be shown (e.g.. Example 3.1 in |KK11| ) that for any e > 

e«(^-^) < Mfc(C), Mfc(C) - X''N{{\})e'^^ < C ^ +oo. 

Moreover, we have, using integration by parts and further applying the Laplace method, 

M,(C) ~ A^Mo(C), C oo, (3.25) 

where 

^o(C)= r (e^''-l-Cv)N{dv). 

J — oo 

Note that the solution ^ (t, x) to (13. 6p satisfies 

e(t,x) = x"'(t)C(xr^^'/'), where ({x) := ^{l,x). (3.26) 

By our choice of A we have Mo(x(t)0 > as t + ^ ^ oo, (t, e S. Therefore, by (K2^ 
and (I3:25|) we have H2: 

M^(t,0 tMo(x(t)0 
Analogously, we have 

^ ^^M2(t,0 ' 

In ((in ((1 V x~\t))M2{t, 0) ) V l) < C(lnlnt + ln(x(t)0) « Int + 
as t + ^ — )■ oo, (t, ^ ^, which provides H3. 
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To show H4, take 6 < and e > such that < fi{[X/f{b), -e]) < oo. Then by flR3|) 
(see Appendix B) we have for \z\ > R with some R large enough 

Q{t,z,R+)>t [ [ (1 - cos{x{t)zuf{s)))fi{du)ds 

J~OQ J{u:0<f{s)u<X} 

> to / (1 — cos{xit)zuf{s)))fi{du)ds 

J-oo Jx/f{b) 

>ctofii[X/fib),-e])\n{xit)e\z\). 

Since we can chose in (1B.3|) c > arbitrary large, condition H4 holds true. Finally, estimate 
(IB.1|) (see Appendix B) provides H5: since /i(R_) > 0, there exists (a, 6) C (— oo,0), q > 0, 
such that < n{[qX / f (b) , X / f (b)]) < oo (note that f{b) < 0), and 



inf e(t,z,[x(t)g,oo))>t inf / / 



(1 — cos{x{t) zuf {s))) fi{du)ds 

/•A//(b) 

> t inf / {{zx{t)uf A l)^ji{du) 
>Cit((x(t)c)2Al). 

Thus, all conditions of Theorem 13.11 are satisfied, and therefore (I3.10p holds true. 

Denote T){x) := H{l,x,^{l,x)). Since ({x) (cf. (K^ ) is the solution to Mi(C(x)) = x, 
we derive from (I3.25P 

({x) ~ Inx, I'(x) ~ —X Inx + A-^x, M2(C(a;)) ~ Ax, x ^ +oo. (3.27) 

From fl3.10p and f l3.27p we deduce that for a given e > there exists y{e) > such that 
( K23\f holds true. □ 



3.4. Properties of Mt{dx) and the completion of the proof. 

Lemma 3.1. For every t > we have Mt{dx) = mt{x)dx with rrit defined by /i3.iy\) . 

Proof. Let g be an arbitrary bounded measurable function. Using the scaling property (13. lip 
of the kernel f{t, s) we can write the definition of the measure Mt in the form 

t r 

H-l/2 . 



g{y)M,{dy) = / / g{t''-'"f{s/t)y)^i{dy)ds 

D J-cx> J{yf{s/t)>X} 

= t / gt{f{s)y)n{dy)ds + t j / gt{f{s)y)^{dy)ds 

J-ooJ{y<-X/f(s)} Jo J{y>X/f{s)} 

:=/_ + /+. 

Here we denote gt{y) '■= g{t^~^^'^y), and recall that /(s) = /(l,s) is monotone on (— oo,0) 
and (0, 1), in particular, / is positive on (0, 1) and negative on (— oo,0). Let us transform 
the integrals J+ and J_ separately. 
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t 



The range of the restriction of / to (0, 1) equals r(H+i/2) ' +°*^ j • Then, making the change 
of variables r = we get 

In the second identity we take into account that, according to f l2.10p . the function i vanishes 
on ^0, r(H+i/2) ) • Then the further change of variables r = y/r and the Fubini theorem give 



Performing similar calculations, we get 



I.=t I gt{r) (I -e\^-] f,{dy) ) dr. 



Adding the expressions for and /_ we get 



/CO POO p - 

g{y)M,{dy) = t / g{t^-^"r)m{r)dr = / g{y) t'/'-^mit^'-'^y) 



dy. 

□ 



Let us summarize: we have fl3.20p and ( I3.23P : in addition, if m is sub-exponential, we 
have fl3.22p by Theorem 3.2 in [K189J. In such a way we obtained the proof of part ii) of 
Theorem 12.21 The following lemma completes the proof of the statement i). 

Lemma 3.2. // yU satisfies fl2.14p . then 



where the constant Ch is defined by f l2.16p . In particular, m G Cd. 
Proof. Write m = m_ + m+. 



m_ r, 



:= r U!i)^^{dy), m+(r):= U(-)^{dy). (3.29) 

J~oo y ^y^ Jo y ^v^ 

On the positive half-axis the function £ can be calculated explicitly: 

Ky) = c^y-(3-^^)/(^-^^) l.>i/r(H+i/2) (3.30) 
with Ch given by fl2.16p . Then 

tn+(r) = c^r-(3"2H)/(i-2H) / y^/(^~^H)^{dy) 



Ch{ I 2/2/(^-'^V(^^y)jr-(3-2^^)/(i-2^^), r^+oo. 

Note that 2/(1 — 2H) > 2 and /i is a Levy measure (that is, J^y^^-^y'^fi{dy) < oo), which 
together with 02.1 4p implies that J^"^ y"^^^^'^^^ fi{dy) < oo. 
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On the negative half-axis one has 

'_,,^-(3-2H)/(l-2If) 



2 r 1--2ff \2/(3-2H) 

with c/f = JZ2H \ 2T{H+i/2) J ' Take arbitrary e > and choose a^, > such that 



Then 

y 



-i(-)Kdy) 



oo 



+ sup (-£{y)) [ (--^(rfy) 



+ / 

r/a^ J~r/bg 

r/h. 

By condition (12. 14^ . one has 
which imphes 

^(3-2H)/(l-2/f) r'\-yf/i^-^")^{dy) -> 0, r ^ +00. 
^ — oo 

Similarly, limr_j.+oo ''"^'^^^~^''^V((~oo, — r]) = and therefore 



r 



r — )■ +00. 



iZ-2H)/il-2H) (]\ < ^(3-2H)/{l-2i/)^^ // _ _^ 

J-r/a, \ yj 

= 6,r2/(^-^^V((-oo,-r/6,]) ^0, 
Then 

limsupr(3-2^)/(i-2^^)m_(r) < (c^ + e) T (-y)'/^'"'''V(^^y)- 

7 >- + 00 J-OO 

The same argument provides the desired lower bound for lim inf^^+oo with ch — s instead 
of Ch + £■ Since e is arbitrary, these two estimates lead to the relation 

° (_^)2/{l-2H)^(^^A -{3-2H)/(l-2H)^ r^+OO, 



which completes the proof. □ 
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4. Two EXAMPLES: THE "EXTREMELY HEAVY- TAILED" CASE 

In this section we give two examples which illustrate the behaviour of pt{x) when condi- 
tion ([2TT1]) fails. In this case we say that the measure fi is "extremely heavy-tailed". 

Example 4.1. Denote by fi-{x) = oo, —x]), fi+{x) = fi{[x, +oo)), x > 0, the "tails" of 
the Levy measure fi, and assume that /i_ and fi^ are regularly varying at +oo, that is, there 
exist a-t G R, and slowly varying functions L±, such that 

see, for example, |Fe71] . Chapter VIII, §8. We investigate the behaviour of the functions m„ 
and m+ introduced in the proof of Lemma 13.21 

Note that a± > 2/(3 — 2H), since otherwise (12. 3p fails. On the other hand, if a±. > 
2/(1 — 2H) then (I2.14p holds true and the required behaviour of m_ and m+ is already 
described in Lemma [3. 2[ Hence, in the sequel we assume that 

«±e f-47T7>^47T7)- (4-1) 



1-2H. 

To simplify the exposition, we exclude from the consideration the critical values a± = 2 / {3 — 
2H) and a± = 2/(1 -2H). 

The asymptotic behaviour of m.+ can be obtained straightforwardly using the standard 
result on the behaviour of the integrals w.r.t. the measures with regularly varying tails, see 
Theorem 2 in |Fe71] . Chapter VIII, §9. The investigation of the behaviour of m.„ is slightly 
more complicated. However, the argument here is quite standard, and therefore we just 
sketch it. 

Write m_ in the form 

m„(r) = ^~'f j{l) f'idy), Hx) := U {^^ . 

It follows from fl3.3ip that there exists a constant C such that $(x) < C(— x)^''^^"^''^^ for 
{—x) small enough and $(x) < C(— x)^/^'^"^'^^ for (— x) large enough. Then, by Theorem 2 
in |Fe71| . Chapter VIII §9, (see also Problem 30 in §10 of the same chapter) we have for A 
small enough and B large enough 

J —oo \T / J— oo 

J-Ar ^"^ J-Ar 

with some explicitly given constant C|^^_ G (0, oo), z = 1, 2. We have 

limsupA2/(i-2^)/^^ = ^2/(1-2^)-- hmsup^^^ = A^/d-^^^)— 
and, similarly, 

hm sup 5^/(3-2^)^^^^ = i?2/(3-2H)-.__ 

r— i>+oo /i_(rj 

Then, by condition (14.11) . for every e > one can choose A and B in such a way that 



14 VICTORIA P. KNOPOVA AND ALEXEY M. KULIK 

Further, the function $ is continuous and positive on [—B,—A]. Therefore, there exists a 
piece-wise constant function such that 

(1 < $ < (1 + (4.3) 

Clearly, one has for every segment (a, b] C R_ 

fi{{ra,rb]) fi_{r) / p_{z)dz, r — )■ oo, with p_(z) = 2;) ^. 

J a 

Therefore, 



f^-[r) J-B 

Combined with (14.21) and (14. 3p . this gives 



^e{^ l^{.dy) ^ j ^e{.z))p^{z)dz, r ^ +00. 



limsup — \— I $ {-\ ddy) <s+ '^^^ / ^iz)p(z)dz. 

One can write in the same fashion the lower bound for lim inf r-jK+oo (we omit the calculation). 
Then, since e > is arbitrary, we finally arrive at 



On the other hand, it follows straightforwardly from Theorem 2 in |Fe71] . Chapter VIII §9, 
that 



m+(r) ~ ^^^i^^ J ^{z)p+{z)dz, r -)■ 00, with p+{z) = a4.z "+ ^ 



Then the function 



1 



r 







m(r) ~ - /i_(r) / ^{z) p^{z) dz + iJ,+ {r) / ^{z) p+{z) dz ] , r +00 



clearly belongs to the class Ld, and thus the statement ii) of Theorem 12.21 holds true. 
Consider, for instance, the "a-stable-like" case 

/i_(r) ~ C_r~", yU+(r) ~ C+r~", with a G ^ 

Then (I2.17P and the above calculations give 



1-2H 



pt[x) 

with 



/ $(z)/i„,c_,c+(c^^), x ^ +00 (4.4) 



fj'a,c-,c+{dz) = a\z\ " ^(^C.1r_{z) + C+1u+{z)^ dz. 

Note that the formal expression for Ha,C-,c+ coincides with that for the Levy measure of an 
a-stable distribution, although for a G (2,2/(1 — 2H)) an "a-stable distribution" itself does 
not exist. 

In contrast to (I2.15p . formula f l4.4p contains explicitly the time parameter t. In addition, 
the polynomial "shape" of the expression in the right-hand side of ( 14. 4 p depends on a, i.e., 
on the "shape" of the tails of the Levy measure fi. 
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Example 4.2. When the Levy measure /i is "extremely heavy-tailed" in the sense explained 
above, the function m may fail to belong to the class JZd at all. Consider the measure 

fc>0 

Then /i(R_) > 0, flCTj) fails, whereas (ES]) is satisfied: 

I \y\'"^'-^''^^l{dy) = 1 + 5^ 22'=/(3-2/f)2fc-2fc/(i-2H) = 1 + J] 2-'- (Wi-f^ < oo. 

k>0 k>0 

Using (13. 30 p . we can write m(r) explicitly: 

n,(^) = _£(_^) + c^,r-(3-2^)/(i-2^^) J2 2^ r>0, 

A;:2fc<rr(_ff+l/2) 

where ch is defined in (I2.16p . To shorten the notation, put c := (T{H + l/2))~^. We have 
for r„ := 2"'c and := (2" — l)c, respectively, 

m(r„) = m(2"c) = -£(-2"c) + CH(2"c)-(3-2H)/(i-2/f) ^ 2^= 

k=0 

and 

n-l 

m(r;) < m((2" - l)c) = -£(-(2" - l)c) + ch((2" - i)c)-(^-m/ii-2H) ^ 2^ 
Note that 

n 

^(-(2" - l)c) ~ £(-2"c), and 2^= ~ 2"+\ as n oo. 

Therefore, 

1 1 m ^ 

m(2"c) - 2' 



fc=0 



fc=0 



and m ^ 

From these two examples one can see that in the "extremely heavy-tailed" case the as- 
ymptotic behaviour of the distribution density of Z/^ is more sensitive with respect to the 
behaviour of the "tails" of the Levy measure fi than in the case where the integrability con- 
dition (I2.14P holds true. If these tails are regularly varying, then (I2.17P holds true with the 
right hand side depending both on t and on the "shape" of the "tails" of fi. On the other 
hand, when the "tails" of fi are both "heavy" and "irregular", the function m may fail to 
belong to the class Ld, which means that we can not state fl2.17p at all. 

Appendix A. Existence of the integral f 1 1.11) 

Proposition A.l. Let < H < 1, H ^ 1/2. Then the integral flLip is well defined for 
every t^Rif, and only if, the Levy measure fi satisfies 112. 3\} . 

Proof. We consider the case < H < 1/2, the calculations in the case 1/2 < H < 1 
are analogous. We check the necessary and sufficient condition for the existence of (11. ip 
formulated in Theorem 2.7 from |RR89) . In our case these conditions can be rewritten as 

f f\s)ds <oo, (A.l) 

J — oo 
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and 



where 
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oo JM. 



(1 A i2{dx)ds < 



oo, 



'2 •- 



(r(/(s)x) - /(s)r(x))Ai(rfx) 



ds < oo, 



TiX 



X, if |x| < 1, 
A, if Ixl > 1, 



(A.2) 
(A.3) 
(A.4) 



and f{s) := /(I, s). Clearly, (lA.ll) is satisfied. We show that a) (\A.2\\ and (12.31) are equiva- 
lent, b) f[0|) follows from flZaj) . 
a) Split 



where 



'11 



'13 



Jl = 111 + /i2 + /l3 + lu, 

»1 

'12 - 



(A.5) 



oo J\uf{s)\<l 

-1 




1 ^l«/(s)l<l 



'14 




1 J \uf (s)\>l 



4t77ttIs|'^ as s — 0—, and 



-oo J\uf{s)\>l 

and estimate the integrals In, i = 1, ..,4, separately. 

Since f{s) 2rctf+i/2) 1^1^'^^^ as s ^ -oo, /(s) - r(//+i/2) 

/(s) = r(H+i/2) (-^ ~ s)'^"^/^ for < s < 1, to check the finiteness of Ii it is enough to substi- 
tute f{s) in the regions (— oo, — 1] and (—1,0) with, respectively, — |s|^~^/^ and — |s|^~^/^, 
and to check the finiteness of the integrals In := /~ Jj^,|<^3/2-H ^12 := /q^ ij^|<,i/2-H 

^13 := /r i|x|>s3/2-« and Ju := /g ij^,|>^i/2-H •••• 
We get: 

r"00 1 



'11 



3-2H 



S 
oo 



r,3-2H 



\x\ fi{dx)ds 

|a;|<s3/2--tf 

\x\^fi{dx)ds + 

\x\<l 



l<|x|<s3/2--ff 



\x\'^fi{dx)ds 



1 



2-2H 
lu- 
lls 



\x\<l 



x\^fi{dx)+ / |a:p/(=^-2^VW 



|x|>l 







1 



X'|<sl/2-H 



x| n{dx)ds < — — / |x| jj,{dx); 



2H ./,,, 



|x|<i 




'14 




1 J|x|>s3/2-tf 

fi{dx)ds 



fi{dx)ds 



(|a;|2/(3-2/f) _ 



X>1 



>/k'|>sl/2-ff 



|^|2/{1-2H)^^^^) ^ / fi{dx). 
\x\<l "'|a;|>l 

Therefore, J < oo if and only if (12. 3 p holds true. 

b) Split h := /21 + /22, where hi := j'^ ii„/(s)|<i - and I22 := ii„/(s)l>i ■ 
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Observe that 

ux — UT{x))fi{dx) = / {ux — - — :)fi{dx)<2 / \ux\fi{dx). 

|a;|<l/|ii| Jl<\x\<l/\u\ fI Jl<\x\<'i-/\u\ 

Then 

J21 < 2 / / \f{s)x\n{dx)ds. (A.6) 

To estimate the right-hand side of (1A.6|) it is enough to estimate 

/21 := r [ -T^f^idx)ds = [ \x\'/^'-^^^^^{dx). 

Thus, (12 .Sp imphes the finiteness of /21, and, consequently, of flA.6p . 
To estimate I22 observe that 

— UT{x))fi{dx) 



|i|>i/|«| F""! 



x|>max(l/|M|,l) F^l FI -''1/|m|<|^|<1 F^ 



"^^ u-^)jj,{dx) + / {- — - — ux)^{dx) 



Then 



/22<2(/ / fi{dx)ds+ [ I \f{s)\fl{dx)ds 

+ f j \f{smdx)ds) 



-1 -'i/l/WI<kl<i 
<Ci{ I I fi{dx)ds+ I \f{s)\ds [ fx{dx) 

1 



+ / / s"^'^'fiidx)ds 

Jo isi/2-H<|a;|<i 

<C2( [ |a:p/(3"2^V(rfa;) + / /i(rfx) + / 

^i|x|>l "'|a;|>l '^|a:|<l 

and the finiteness of the right-hand side is imphed by (I2.3p . □ 



Appendix B. Existence of the distribution density 

Proposition B.l. Let < H < 1, H ^ 1/2. Then, under and our standing assumption 
Iji2.4\ ), the integral U.l\) possesses for any t ^ the distribution density pt G C^(R). 

Remark B.l. In the non-Markov case H ^ 1/2, the kernel /(t, s) provides a strong "smooth- 
ifying" effect in the sense that the weakest possible non- degeneracy assumption (12. 4p is al- 
ready sufficient for the integral (11 .ip to possess a smooth distribution density. We refer to 
[KKllj . Section 3, for the detailed discussion of various forms of the "smoothifying" effect 
for Levy driven stochastic integrals with deterministic kernels. 



For the proof we use the following statement, see Lemma 3.3 in |KK11| . 
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Proposition B.2. a) For a positive function h{s) having a continuous non-zero derivative 
on some interval [a, b] C R., one has 

j {l-cos{h{s)x))ds>c{x^M). (B.l) 

J a 

b) For a positive convex on (—00,6) C R function h{s), satisfying 

lim e^^'h{s) = +00 for all 7 > 0, (B.2) 

one has ^ 

/ {l-cos{xh{s)))ds>c\n\x\ (B.3) 
J —00 

for all c > and \x\ big enough. 

Proof of Proposition \B. 1[ Recall (cf. ( 12. 5p ) that the characteristic function of Zf^ is of the 
form (j){t,z) = e'^^^'~^K For a fixed t, the function h{s) = —t^~^^'^f{s) satisfies (1B.2P with 
6 = 0. Since /i(R) > 0, there exists g > such that 

Q := max{/i((-oo, -g]), 00)} > 0. 

Then using flB.3|) for \z\ large enough we get 

-Re ^(t, -z)>t f I (1 - co&{t^-^''^ f{s)uz))^i{du)ds > tcQ In \qt"-^/h\. 

J — 00 J \u\>q 

Since c > is arbitrary, the function |</)(t, = e"R'e'i'(i,-z) jg jntegrable in z for any a > 0. 
Therefore the density pt is well defined and belongs to as the inverse Fourier transform 
of z): 

Pt{^) = ^ I e-'^^(l>{t,z)dz. (B.4) 

□ 



Remark B.2. Literally the same argument implies that, if the truncation level A > is 
chosen small enough, then = e^^'^'^i^*'^) is integrable in z for any a > 0, t > 0, see 

f l3.12p . This implies the existence ofpt{x), see f l3.15p . 
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